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Page 5, first paragraph under "DISCUSSION AND RESULTS": The sentence that begins 
with "Horizontal wind and vertical wind profiles..." should read 
"Horizontal wind and vertical wind profiles in figure 7 were taken at h = 
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SUMMARY 


A simple downburst model has been developed for use in batch and real- 
time piloted simulation studies of guidance strategies for terminal area 
transport aircraft operations in wind shear conditions. The model 
represents an axisymmetric stagnation point flow, based on velocity profiles 
from the Terminal Area Simulation System (TASS) model developed by Proctor 
[ref. 3,4] and satisfies the mass continuity equation in cylindrical 
coordinates. Altitude dependence, including boundary layer effects near the 
ground, closely matches real-world measurements, as do the increase, peak, 
and decay of outflow and downflow with increasing distance from the 
downburst center. Equations for horizontal and vertical winds were derived, 
and found to be infinitely differentiable, with no singular points existent 
in the flow field. In addition, a simple relationship exists among the 
ratio of maximum horizontal to vertical velocities, the downdraft radius, 
depth of outflow, and altitude of maximum outflow. In use, a microburst can 
be modeled by specifying four characteristic parameters. Velocity 
components in the x, y, and z directions, and the corresponding nine partial 
derivatives are obtained easily from the velocity equations. 


INTRODUCTION 

Terminal area operation of transport aircraft in a windshear 
environment has been recognized as a serious problem. Studies of aircraft 
trajectories through downbursts show that specific guidance strategies are 
needed for aircraft to survive inadvertant downburst encounters. In order 
for guidance strategies to perform in simulations as in actual encounters, a 
realistic set of conditions must be present during development of the 
strategies. Thus, airplane and wind models that closely simulate real-world 
conditions are essential in obtaining useful information from the studies. 

Wind models for use on personal computers, or for simulators with 
limited memory space availability, have been difficult to obtain because 
variability of downburst characteristics makes analytical models 
unrealistic, and large memory requirements make use of numerical models 
impossible on any except very large capacity computers. 

Bray [ref. 1] developed a method for analytic modeling of windshear 
conditions in flight simulators, and applied his method in modeling a 
multiple downburst scenario from Joint Airport Weather Studies (JAWS) data. 
However, the altitude dependence of his model is not consistent with 
observed data, and, although flexibility in sizing the downbursts is built 
into the model, it does not maintain the physical relationships which are 
seen in real-world data among the sizing parameters. In particular, 
boundary layer effects should cause radial velocity to decay vertically to 
zero at the ground, as does the vertical velocity. 

In a study conducted at NASA Langley Research Center, three different 
guidance strategies for a Boeing 737-100 airplane encountering a microburst 
on takeoff were developed [ref. 2] . These strategies were first developed 
using a personal computer, and then implemented in a pilot-in-the-loop 
simulation using a very simple wind model in both efforts [fig. 1] . This 
model consisted of a constant outflow outside of the downburst radius and a 
constant slope headwind to tailwind shear across the diameter of the 
downburst. It was recognized that a more realistic wind model could 
significantly alter the outcome of the trajectory. For the subsequent part 



of this study, which involves altering the airplane model to simulate 
approach to landing and escape maneuvers and additional takeoff cases, a 
more realistic wind model was preferred. The simple analytical model 
outlined in this report was developed for this purpose. 
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SYMBOLS 

Joint Airport Weather Studies 

Northern Illinois Meteorological Research On Downbursts 
radius of downburst shaft (ft) 

radial coordinate (distance from downburst center) (ft) 

Terminal Area Simulation System 

velocity in r-direction (or x-direction) (kts) 

velocity in y-direction (kts) 

velocity in z-direction (kts) 

magnitude of maximum vertical velocity (kts) 

magnitude of maximum horizontal velocity (kts) 

horizontal (runway) distance, airplane to downburst center (ft) 

horizontal (side) distance, airplane to downburst center (ft) 

airplane altitude above ground level (ft) 

depth of outflow (ft) 

height of maximum U-velocity (ft) 

height of half maximum U-velocity (ft) 

characteristic height, out of boundary layer (ft) 

characteristic height, in boundary layer (ft) 

scaling factor (s” 1 ) 


DEVELOPMENT OF VELOCITY EQUATIONS 

Beginning with the full set of Euler and mass continuity equations, 
some simplifying assumptions about the downburst flow conditions were made. 
Effects of viscosity were parameterized explicitly, and the flow was assumed 
to be invariant with time. The downburst is axisymmetric in cylindrical 
coordinates, and characterized by a stagnation point at the ground along the 
axis of the downflow column. The flow is incompressible, with no external 
forces or moments acting on it. 
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The resulting mass conservation equation is 


V- V= 0 . 


( 1 ) 


Written out in full, equation 2 is 


ilL+iw + U..Q 
3r dz r 


This equation is satisfied by solutions of the form 


w = g (r )q(z) 


f(r 2 ) 


u = — ~ — p(z) 


provided that 

f'(r 2 ) =lg(r 2 ) 
q'(z) = - \p(z) 4 


( 2 ) 


(3a) 

(3b) 


(4a) 

(4b) 


f '( r 2 ) = 

Note that dr . To solve this system of equations, solutions were 

assumed for two of the functions and the other two were obtained from 
equations 4a and 4b. 

It was desired that the velocity profiles of this analytic model 
exhibit the altitude and radial dependence shown in the large-scale 
numerical weather model TASS (Terminal Area Simulation System) [ref. 3,4]. 
The TASS model is based on data from the Joint Airport Weather Studies 
(JAWS) [ref. 5], and provides a three-dimensional velocity field, frozen in 
time, for given locations of an airplane within the shear [ref. 6] . Figure 
2 shows dimensionless vertical profiles of horizontal velocity, u, for TASS 
data, laboratory data obtained by impingement of a jet on a flat plate, and 
data from NIMROD (Northern Illinois Meteorological Research On Downbursts) 
[ref. 7] Specific points of interest are the maximum horizontal velocity 
(located 100 - 200 meters above the ground), below which is a decay region 
due to boundary layer effects, zero velocity at the stagnation point on the 
ground, and an exponential decay with altitude above the maximum velocity 
altitude. Vertical velocity profiles from TASS data are shown in figure 3, 
also exhibiting a decay to zero at the stagnation point. 

The radially varying characteristics desired for the horizontal wind 
were two peaks of equal magnitude and opposite direction located at a given 
radius, with a smooth, nearly linear transition between the two. Beyond the 
peaks, the velocity should show an exponential decay to zero. The vertical 
velocity was required to have a peak along the axis of symmetry (r=0) , and 
decay exponentially at increasing radius. 

A pair of shaping functions that gave velocity profiles matching TASS 


3 



data as required are given below. 


? (r 2 ) = 


HrM) 


p( z) = e _z/z * - e” z/£ 


The remaining solutions were found by integrating equations 4a and 4b, 
yielding : 
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1 - e 


-(r/R) 


fir 2 ) 

q( z) = — X.{ e(e -z/E - ]) - z*(e~ z/z * - l)} 


Figures 4 and 5 show plots of these shaping functions. 

Combining the functions as in equation 3, the horizontal and vertical 
velocities are expressed as 


u = 


Xr 

2 r 


w = - Xe" (r/R)2 [£(e- z/£ - l) - z*(e _z/z * - l)] . 


(5) 


( 6 ) 


By taking derivatives of equations 5 and 6 with respect to r and z, 
respectively, and substituting in equation 2, it can be shown that the 
velocity distributions satisfy continuity. 

The parameters z* and e were defined as characteristic scale lengths 
associated with "out of boundary layer" and "in boundary layer” behavior, 
respectively. Analysis of TASS data indicated that z* = z m2 , the altitude 
at which the magnitude of the horizontal velocity is half the maximum value. 

It was also noted that the ratio 



0.22 


To determine the location of the maximum horizontal velocity, the 
partial derivatives of u with respect to r and z were set equal to zero. 
The resulting equation for the r-derivative is 

<*)’-• w ->. 

The resulting equation for the z-derivat ive is 



1 

(z*/e) - 1 


ln(z * /e) 


4 



Recalling that z m /z* = 0.22, the values 1.1212 and 12.5 were obtained from 
iteration for the ratios r/R and z*/e, respectively. 

Using these values, the maximum horizontal velocity can be expressed 
as u max - 0.2357 X R. The maximum vertical wind is located at r - 0 and z = 

z h , by definition, and is given by w max = Xz* (e“ < z h /z *> -0.92) . 

A ratio of maximum outflow and downflow velocities can be formed 


0.2357R 

-(Z /z*) 

z* (e - 0.92) 


The scaling factor, X, was determined by using either of equations 5 
or 6 for horizontal or vertical velocity, and setting it equal to the 
maximum velocity, u max or w max , respectively. Solving for X resulted in: 


■—( z / z* ) 

z* (e - 0.92) 0.2357R 

The velocity equations were easily converted to rectangular 
coordinates, as shown in the Appendix. Partial derivatives with respect to 
x, y, and z were obtained by differentiating the velocity equations, and are 
also listed in the Appendix. 


DISCUSSION AND RESULTS 

Vertical and horizontal velocity profiles for u and w are shown in 
figures 6 and 7. Four profiles are shown for each component. The 
horizontal wind profiles in figure 6 were taken at the radius of peak 
outflow (r = 1.1212 R) and at about one-fourth that radius (r = 0.3 R) , 
where the maximum outflow is approximately half the value at the peak 
outflow radius. The vertical wind profiles were taken at the radius of peak 
downflow (r = 0) and at r = 0.3 R. Horizontal wind and vertical wind 
profiles in figure 7 were each taken at two altitudes, h = z m (altitude of 
maximum outflow) and h = z h (depth of outflow) . The outflow velocity is 
reduced by half (from 37 kt to 18 kt) , while the downflow velocity is 
increased by about 38 percent (from 21 kt to 29 kt) in the same distance. 

This analytical model is compared with TASS, laboratory, and NIMROD 
data in figure 8. The figure shows that, when nondimensionalized by the 
altitude of half-maximum outflow (z*) and by the maximum outflow (u = u max ) , 
the analytical model agrees closely with the other data. 

Different shears can be modeled by specifying four parameters, and the 
location of downburst center relative to the airplane flying through it. 

The four parameters are: 1) a characteristic horizontal dimension; 2)maximum 
wind velocity; 3) altitude of maximum outflow; and 4) depth of outflow. The 
characteristic horizontal dimension specified is the radius of the downdraft 
column, noting that this is about 89 percent of the radius of peak outflow. 
The maximum wind velocity can be either horizontal or vertical. 
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CONCLUDING REMARKS 


The analytic microburst model developed for use in real-time and batch 
simulation studies was shown to agree well with real-world measurements for 
the cases studied. The functions chosen for the model showed boundary-layer 
effects near the ground, as well as the peak and decay of outflow at 
increasing altitudes, and increasing downflow with altitude. The 
exponential increase and decay of downflow and outflow (in the radial 
direction) are also characterized by the model. Equations for horizontal 
and vertical winds are simple and continuously differentiable, and partial 
derivatives in rectangular or cylindrical coordinates can be easily obtained 
by direct differentiation of the velocity equations. The governing equation 
for this system is the mass conservation law, and the analytic velocity 
functions developed here satisfied this condition. The model is sustained 
by a strong physical basis and yields high fidelity results, within the 
limitations of maintaining simplicity in the model, and variability of the 
microburst phenomenon. Parameterization of some of the characteristic 
dimensions allows flexibility in selecting the size and intensity of the 
microburst . 
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APPENDIX 


Define intermediate variables to simplify written equations: 
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From TASS 


Other Relationships 


■= 0.22 


= 12.5 


Maximums 


W v = 0.2357X.R 

max 

w =w x 

1 max max 

-<z /z*) 

W = Xz*e 

h 


- 0.92) 


( X is determined from the above relationships) 
w 

X max 0.2357 R 


W ” Z K /Z * 

h max z* (e -0.92) 


Variable List 


w w x 

z* * altitude where x is half the value of max (ft) 

£ *= characteristic height of boundary layer effects (ft) 
z 

h = depth of outflow (ft) 

2 

m « altitude of maximum outflow (ft) 

. -l 

A. scaling parameter ( 3 ) 

r * radial distance from airplane to downburst (ft) 
h * altitude of airplane (ft) 

R - radius of downdraft (ft) 

X y 

ad ad = x,y coordinates, airplane to microburst (ft) 
W x 9 W y ' W h 

max max max maximum winds, x, y, and h directions 


i 
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Figure 2 Vertical Profile of Microburst 
(Nondimensional) 
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Figure 3 Vertical Profile of Microburst Downflow 
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Figure 7 Radial Velocity Profiles For Analytical Model 
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Figure 8 Comparison of Wind Model Vertical Profiles 
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